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Abstract 



| Random sets with long-range dependence can be generated using a Boolean 

model with power-law grain sizes. We study thinnings of such Boolean mod- 
els which have the hard-core property that no grains overlap in the resulting 
germ-grain model. A fundamental question is whether long-range dependence 
is preserved under such thinnings. To answer this question we study four nat- 
ural thinnings of a Poisson germ-grain model where the grains are spheres 
with a regularly varying size distribution. We show that a thinning which 
favors large grains preserves the slow correlation decay of the original model, 
whereas a thinning which favors small grains does not. Our most interesting 
finding concerns the case where only disjoint grains are retained, which corre- 
sponds to the well-known Matern type I thinning. In the resulting germ-grain 
00 ■ model, typical grains have exponentially small sizes, but rather surprisingly, 

the long-range dependence property is still present. As a byproduct, we ob- 
s ! , tain new mechanisms for generating homogeneous and isotropic random point 

configurations having a power-law correlation decay. 

o 

1 Introduction 

Consider a random closed set which can be expressed as a union of compact sets 
in the d-dimensional Euclidean space M. d . The compact building blocks of the 
random set are called grains, the collection of grains germ-grain model, and the 
union of grains grain cover. A germ-grain model is called hard-core if the grains are 
disjoint with probability one. Hard-core germ-grain models (a.k.a. random packing 
models) provide an important class of mathematical tools for the natural sciences, 
allowing to model and analyze the statistical features of disordered porous materials 
[OMOO, SSW02]. Besides natural sciences, these models have found applications in 
engineering when analyzing the performance of medium access protocols in wireless 
data networks (e.g. [BB09, Haell, NB12]). 
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A key statistical feature of a random set is its covariance function, which de- 
scribes how much more or less likely it is to find matter at a given distance from 
a location containing matter, compared to finding matter in an arbitrary location. 
While most germ-grain models studied in the literature have a rapidly decaying 
covariance function, certain experimental studies in astronomy [JMST05] and ma- 
terials science [SSW02] display real-world data where the statistically estimated 
covariance function appears to decay exceptionally slowly, following a power law 
with some exponent j5 > for large distances r. When (3 < d, such models are 
long-range dependent in the sense that 

r Var(|Xni? r |) 

hm sup 2 = oo, (1.1) 

where \X n B r \ denotes the volume of the region covered by the random set X 
within the closed ball B r with radius r centered at the origin [DVJ08, Sec 12.7]. 
Long-range dependence causes anomalous behavior to several statistical features of 
the model, as is well understood in time series analysis [Sam06]. Note that for a 
homogeneous random set in dimension d = 1, property (1.1) is equivalent to the 
usual notion of long-range dependence, 

limsup = oo, 

n— >oo Tl 

of the time series X^ = \X H {k — 1, k]\. 

Our goal in this article is to construct parsimonious germ-grain models having 
the hard-core and long-range dependence property. In the presence of long-range 
dependence, the requirement of parsimony, i.e. having a small number of model 
parameters, is especially important because long-range dependence tends to reduce 
the robustness of the statistical estimators of model parameters [CSN09]. Long- 
range dependent germ-grain models are easy to generate using a Boolean model — a 
germ-grain model with random power-law distributed sizes and independently and 
uniformly scattered centers — but the resulting model is not hard-core by construc- 
tion. To make it hard-core, we shall follow Matern's approach [Mat60] of thinning 
out a selected collection of overlapping grains from the proposed Boolean model 
so that the resulting collection of grains is disjoint. Whether this approach is fea- 
sible for obtaining hard-core models with long-range dependence depends on the 
following question: 

Is the power-law covariance decay of the proposed Boolean model pre- 
served after making it disjoint by thinning? 

To answer the above question, we shall analyze in detail the following natural 
thinning mechanisms: 

• Large retained. Let the thinned model consist of those grains in the original 
Boolean model which are not overlapped by any larger grain in the original 
model. 
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• Random retained. Assign independent random weights to the grains. Let 
the thinned model consist of those grains in the original model which are 
not overlapped by any heavier grain in the original model. (This thinning 
corresponds to Matern type II.) 

• Small retained. Let the thinned model consist of those grains in the original 
model which are not overlapped by any smaller grain in the original model. 

• Isolated retained. Let the thinned model be the set of grains in the original 
model which do not overlap with any other grain in the original model. (This 
thinning corresponds to Matern type I.) 

We remark that — unlike the Matern type III hard-core model [NB12] — the above 
thinnings are local in that the decision whether a proposed grain shall be retained 
or not is made solely by looking at the grains which intersect it. 

For simplicity, we shall restrict to spherical models where the grains are closed 
balls. Figure 1 illustrates the above four thinnings applied to a simulated sample 
of a Boolean model in M 2 where the grain centers have mean density A = 0.05 and 
the grain radii have a Pareto distribution F(r) = 1 — r~ a , r > 1, with tail exponent 
a = 2.5. 



Original Large retained Random retained 




-40 -20 20 40 -40 -20 20 40 



Figure 1: Original model and its four hard-core thinnings, where the original model 
is generated by disks having a Pareto distribution with tail exponent a = 2.5. 

The above thinnings will be analyzed collectively by viewing them as instances 
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of a general weight-based thinning mechanism, following Mansson and Rudemo 
[MR02]. Standard formulas of Palm calculus allow to write down closed-form ana- 
lytical formulas for the radius distribution of a typical grain, the covariance function 
of the grain cover, and the two-point correlation function of the grain centers for 
general hard-core germ-grain configurations generated by weight-based thinnings. 
Using the theory of regular variation, we analyze the long-range behavior of these 
quantities under the assumption that the grain radii in the proposed Boolean model 
follow a power-law distribution with tail exponent a > d. 

The main results of this article (Theorems 5.1, 6.1, 7.1, 8.1) are summarized in 
Table 1 below. From the table, we can draw the following conclusions: 



Model 


Radius distribution 
of a typical grain 


Covariance function 
of grain cover 


Correlation function 
of grain centers 


Original 
Large retained 
Random retained 
Small retained 
Isolated retained 


power law (a) 
power law (a) 
power law (a + d) 
exponential 
exponential 


power law (a — d) 
power law (a — d) 
power law (a — d) 
exponential 
power law (a — d) 


zero 

power law (a — d) 
power law (a — d) 
exponential 
power law (a — d) 



Table 1: Long-range decay of key statistical characteristics of the original Boolean 
model and the hard-core germ-grain models obtained by thinning. 



• The power-law covariance decay and long-range dependence (when a < 2d) 
of the grain cover are preserved under all thinnings except small retained. 

• Whereas the random point configuration (a.k.a. point process) of grain centers 
in the proposed Boolean model is completely uncorrelated, the corresponding 
point configurations in all thinned models except small retained have a power- 
law two-point correlation function. 

• The heavy tail of the grain radius distribution is destroyed by small retained 
and isolated retained thinnings. The other two thinnings preserve the power- 
law structure of the tail distribution: under large retained with the same 
exponent, under random retained with a larger exponent corresponding to a 
lighter tail. 

Table 1 also reveals a striking feature of the isolated retained thinning mechanism: 
The resulting grain cover and the resulting point configuration of grain centers both 
exhibit long-range dependence although the grain size distribution is light-tailed. 
This seemingly paradoxical phenomenon can be explained by inspecting the empty 
space: Any region of space not covered by the thinned germ-grain model is likely 
to have been contained in a big grain of the proposed model that was removed in 
the thinning, and therefore, a large neighborhood of this empty region is likely to 
be empty, too. 

This article may be seen as a continuation of the works of Mansson and Rudemo 
[MR02] and Andersson, Haggstrom and Mansson [AHM06], who analyzed first- 
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order statistical properties of hard-core germ-grain models obtained by weight- 
based thinnings. In [MR02, Cor 3.1] it was also shown that large retained thinning 
preserves the tail behavior of the typical grain radius whenever the proposed grain 
radius distribution is continuous. A slightly more general thinning framework was 
recently introduced by Nguyen and Baccelli [NB12], who derived differential equa- 
tions characterizing the generating functional of the random point configuration 
formed by the thinned grain centers. Earlier work on the covariance analysis of 
random sets includes Bohm and Schmidt [BS03], who derived a short-range ap- 
proximation for the covariance function of a general homogeneous random set. 
Snethlage, Martinez, Stoyan, Saar [MVDE02] (see also references therein) pro- 
vide a nice summary of random point configuration models where the two-point 
correlation function has a power-law behavior on short distances. Earlier works 
on long-range dependent random sets appear mostly restricted to random point 
configuration in dimension d = 1. Among these, Daley and Vesilo [DV97] estab- 
lished the following elegant preservation property for many queueing systems: the 
point configuration of the departure times is long-range dependent if and only if 
the same is true for the arrival times. Daley [Dal99] showed that a renewal point 
process is long-range dependent if the interpoint distances have an infinite second 
moment, and Kulik and Szekli [KS01] extended this observation to one-dimensional 
point configurations with positively associated interpoint distances. Vamvakos and 
Anantharam [VA98] showed that the long-range dependence of a point process is 
preserved by a leaky bucket flow control mechanism for data traffic. A study focused 
on the long-range dependence of multidimensional random sets is the recent work 
of Demichel, Estrade, Kratz, and Samorodnitsky [DEKS11], who studied whether 
random sets having power-law decaying chord length distributions, closely related 
to the covariance function of the random set, can be generated as a level set of a 
Gaussian random field — they found that in wide generality (merely assuming that 
the underlying Gaussian field is mixing), this is not possible. 

Let us summarize the notational conventions used in this paper. The symbol P 
stands for the probability measure on some abstract probability space which governs 
all randomness in the models, and E,Var,Cov denote the expectation, variance, 
and covariance with respect to P, respectively. The symbol B r (x) denotes the 
closed unit ball with center x and radius r in the <i-dimensional Euclidean space 
Mr. We use B r as shorthand for B r (o), where o is the origin of R d . For a Borel 
set B in we denote by \B\ its Lebesgue measure, and by 1b(x) or l(x € B) 
its indicator function. The symbols dx, dy, etc. refer to the Lebesgue measure 
in M. d . The symbol M+ denotes the positive real numbers including zero. The 
symbol F(dr) refers to integration with respect to a probability measure F on 
R_l_, whereas F(r) = F[0,r] and F(r) = 1 — F(r) stand for the corresponding 
cumulative distribution function and the complementary cumulative distribution 
function, respectively. The minimum and maximum of real numbers a and b are 
denoted by a A b and a V b, respectively. When convenient, we denote J a °° = Jj a s, 

fo = J[o 61' an d fa = f(ab]^ OI ® <a< b<oo. For functions / and g defined on the 
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positive real line, we denote / ~ g if f(t)/g(t) — > 1 as i — > oo. 

The rest of the paper is organized as follows. Section 2 summarizes prelimi- 
naries on random Boolean models needed later in the text. Section 3 introduces 
a weight-based thinning mechanism which produces hard-core germ-grain models 
from Boolean models and list formulas for the second-order statistics of the models 
so obtained. Section 4 contains a long-range analysis of the second-order statistics 
of the previous section. The main results of Table 1 are proved case- by-case in Sec- 
tion 5 (isolated retained), Section 6 (random retained), Section 7 (large retained), 
and Section 8 (small retained. Section 9 concludes the paper. 

2 Boolean models with power-law grain radii 

A spherical Boolean model is a random collection of closed spheres, where the sphere 
centers are independently and uniformly scattered in M. d and the sphere radii are 
independent and identically distributed random variables in M+. Mathematically, 
a spherical Boolean model can be defined as a Poisson random measure <3? on W d x 
M + with intensity measure XdxF(dr), where A is a positive constant and F is 
a probability measure on M + such that f r d F(dr) < oo. We identify each pair 
(x,r) € <3? with the closed ball B r {x) with center x and radius r and — conforming 
to the terminology of more general germ-grain models — such pairs wills be called 
grains. The random closed set 

X= |J B r (x) 

is called the grain cover of <3?, and we denote by 

& g = {x eR d : (x, r) E <3? for some r} 

the random point configuration in M. d formed by the grain centers of 3>. Note that 
$ g is a homogeneous Poisson random measure on M rf with intensity measure Xdx. 
The parameter A thus equals the mean density of grain centers, and the probability 
measure F is the common distribution of grain radii. For general definitions and 
details about random sets and random measures, see for example [DVJ08, Mol05, 
SW08, SKM95]. 

The covariances of the random set X are denoted by k(x, y) = Cov(lx(x), lx(y)), 
where lx is the indicator function of X. Because the distribution of X is shift- 
invariant by construction, the covariances are given by k(x,y) = k(x — y), where 
the covariance function k(z) = k(o,z) is given by the well-known formula (e.g. 
[SKM95, Sec. 3.1]) 

k(z) = (1 - pf (e A / \Mo)nB r{ z)\ F(dr) _ ^ (21) 

and where p is the volume fraction of X given by 

p =l- e -*S\Br\F{dr)_ 
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Formula (2.1) indeed shows that k{z) depends on z only through \z\, which is 
evident because X is isotropic by construction. Using this formula we may also 
deduce that 

k(z) ~ A(l-p) 2 J \B r (o)C\B r (z)\F(dr) as \z\ -> oo, 

where we denote f(z) ~ 5(2) if f(z)/g(z) — > 1 as |z| — > 00. When the grain 
radius distribution F follows a power law with tail exponent a > d, so that F(r) = 
1 — £(r)r~ a for some slowly varying function £ (see Appendix A for details), it 
follows by using Lemma B.l in Appendix B that 

k(z) ~ A(l - pfc a4 £(\z\)\z\-^- d ^ as |z| -> 00. 

Thus, when the radius distribution follows a power law with tail exponent a > d, 
then the covariance function k(z) follows a power law with tail exponent a — d. 
Especially, the Boolean grain cover X is long-range dependent in the sense of (1.1) 
for a £ (d, 2d). 

3 Weight-based thinning 

In this section we shall study a weight-based thinning mechanism which maps a 
Boolean model into a hard-core germ-grain model consisting of nonoverlapping 
grains [MR02, NB12]. This thinning mechanism is defined by assigning random 
weights to the grains of the Boolean model, and retaining those grains which are 
not overlapped by any other grain in the Boolean model with a higher or equal 
weight. 

3.1 Thinning mechanism 

A weighted spherical Boolean model is defined as a Poisson random measure $ on 
M. d x R + x M + with intensity measure 

A(dx,dr,dw) = XdxF(dr)G r (dw), 

where A > 0, F is a probability measure on R + such that f r d F(dr) < 00, and G is 
a probability kernel on R + (a family of probability measures G r on M + indexed by r 
such that r 1— > G r (A) is measurable for measurable A C A triplet (x, r,w) € $ 
is identified as a grain with center x, radius r, and weight w. As in Section 2, the 
constant A is the mean density of grain centers and the probability measure F is the 
distribution of grain radii. The probability measure G r is the weight distribution 
of a grain with radius r. 

We say that two distinct grains are neighbors if they intersect each other, and 
we denote the set of neighbors of a reference grain (x, r, w) by 

N x>r ,w = {(a/.r'X) G R d x R + x R + \ {(x,r,w)} : B r /(x') fl B r {x) + 0} . (3.1) 
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The thinning of a weighted spherical Boolean model $ is now defined by = 
T(3>), where 

T($) = {(x,r,w) G $ : io > w' for all (x',r',u/) G $ n N x ^ w } . (3.2) 

To rephrase the definition, we say that a grain (x',r',w') obstructs grain (x,r,w) 
if (x',r',w') is a neighbor of (x,r, w) and w' > to. Then by definition, the thinned 
germ-grain configuration $ th consists of grains in $ which are not obstructed by 
any other grain in Note that two overlapping grains with equal weights obstruct 
each other, and will be both removed. 

The following choices of G r yield the four thinnings which shall be analyzed in 
detail in Sections 5-8. 

• Large retained. The weight of each grain is set equal to its radius, so that 
G r {alw) = 5 r (dw). 

• Random retained. The grains are assigned independent uniformly distributed 
random weights, so that G r (dw) = lr 0)1 \(w)dw. 

• Small retained. The weight of each grain is set equal to the inverse of its 
radius, so that G r (dw) = 5i/ r (dw). 

• Isolated retained. All grains are assigned weight one, so that G r (dw) = 
Si(dw). 

3.2 Retention probability 

The retention probability of a reference grain (x, r, w) is defined as the probabil- 
ity that (x,r, w) belongs to the germ-grain configuration obtained by thinning the 
union $ U {(x,r,w)}. Because this probability does not depend on x (see Proposi- 
tion 3.1), we shall denote it by 



The quantity h(r, w) may be regarded as the probability that a typical grain with 
radius r and weight w in the proposed Boolean model is retained (see e.g. [SW08, 
SKM95]). Analogously, the weight-averaged retention probability 



may be regarded as the probability that a typical grain of radius r in the proposed 
Boolean model is retained. The following result [MR02, Thm. 2.2] gives a formula 
for the retention probability. For the reader's convenience we will include the proof 
here. 



h(r, w) = P((x, r, w) <E T($ U {(x, r, w)})). 




(3.3) 
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Proposition 3.1. The retention probability of an arbitrary reference grain (x, r, w) S 
R rf x M + x M + does not depend on x, and is given by 

h(r,w) = expj-A J \B r+s (o) \ G s [w, oo) F(ds)j . (3.4) 

Proof. Fix a reference grain (x, r, w) and denote <3?' = <3? U {(x, r, w)}. By definition, 
the reference grain belongs to the thinned configuration T(<3?') if and only if w > w' 
for all (x',r',w') G N x r ,v) n where iV^. r w is the neighbor set of (x, r, w) defined 
by (3.1). Observe that N X)r>w r\$' = N Xt r iW n<&, because no grain is its own neighbor 
by definition. As a consequence, the retention probability can be expressed using 
the the intensity measure of the Poisson point configuration $ according to 

P((x,r,w) e T($')) = PmA X;r!W ) = 0) = e - A (^.-»), 

where 

Ao,r,w = {{x\r',w') G N X) r lW : w' > w} 
is the set of grains obstructing (x,r,w). The claim now follows because 

A(A X ^ W )= f [ [ l(\x-x'\<r + r')l(w'>w)G r/ (dw')F(dr')Xdx' 

JR d JR+ Jr+ 

= X \B r+r ,{o)\G r ,[w,oo)F(dr'). 
Jr + 

□ 



3.3 First-order statistics of the thinned model 

Let us summarize some key formulas about the first-order statistics of the thinned 
germ-grain model <3? th which were obtained in [AHM06, MR02]. The mean density 
of grain centers in the thinned model is given by 



Ath = A/ h(r)F(dr), (3.5) 

where h(r) is the weight-averaged retention probability defined in (3.3), and the 
radius distribution of a typical grain in the thinned model equals 

A 1'°° 

F th (r) = l-— / h(s)F(ds). (3.6) 

Moreover, the volume fraction of the thinned grain cover 

X th = |J B r (x) 
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is given by 

Pth = \[ \B r \h{r)F(dr). (3.7) 

Note that the quantity J h(r)F(dr) in (3.5) may be regarded as the probability 
that a randomly chosen grain in the proposed Boolean model is retained by the 
thinning mechanism. 

3.4 Pair retention probability 

The pair retention probability of a given pair of reference grains {xi,r\,w\) and 
(%2, r 2, w 2 ) is defined as the probability that both reference grains belong to the 
germ-grain configuration obtained by thinning the union $' = <J>U{(xi, r±, wi), (x 2 ,r 2 ,w 2 )}- 
Because this probability depends on x\ and X2 only through their distance (see 
Proposition 3.2), we shall denote it by 

h 2 (u,n,wi,r 2 ,w 2 ) = P({(xi,ri,u;i), (x 2 , r 2 , w 2 )} € T($'))> (3.8) 

where u = \x\ — x 2 \. The weight-averaged pair retention probability is defined by 

h 2 (u,r 1 ,r 2 )= / h 2 (u,ri,w 1 ,r 2 ,w 2 )G n (dwi)G r2 (dw 2 ). (3.9) 

Proposition 3.2. The pair retention probability of two reference grains (xi,ri,w{) 
and (x 2 ,r 2 ,w 2 ) depends on x\ and x 2 only through the distance u = \x\ — x 2 \. For 
u < r\ + r 2 this probability equals zero, and for u > n + r2, 

h 2 (u, ri , Wl ,r 2 ,w 2 ) = h(r 1 ,w 1 )h(r 2 ,W2)e T ( u ' r ^ r ^ 

where h{r\,W\) and h(r 2 ,w 2 ) are the retention probabilities defined by (3.4), and 

T(u,ri,wi,r 2 ,w 2 ) = X / \B r+ri (xi) n B r+r2 (x 2 )\ G r [w\ V w 2 , 00) F(dr) 

JR + 

is the mean number of grains in <I> which simultaneously obstruct both reference 
grains. 

Proof. Fix two reference grains (aci, r±, Wi) and (x 2 ,r 2 ,w 2 ) and assume that they do 
not overlap, so that |xi — a^l > ^1 + ^2- Denote = $ U {(xi, n, wi), (x2, r2, ^2)}- 
Recall that grain (xi, ri, u>x) belongs to T($') if and only if w\ > w for all (x, r, w) € 
N(xi, r\, wi) PI Because no grain is its own neighbor by definition, and be- 
cause the two reference grains are not neighbors, we see that N(xi, r±, w\) n = 
N(x\,r\,w{) H By symmetry, a similar conclusion also holds for the other refer- 
ence grain. 

We conclude that for i = 1,2, grain (xj, ri,Wi) is retained if and only if <3?(Aj) = 
0, where 

Ai = {(x,r,w) E N(xi,n,Wi) : w >Wi} 
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is the set of grains obstructing (xj, r^, Now the pair retention probability can 
be written as 

h 2 = P($(A 1 UA 2 ) = 0). (3.10) 
The number of grains in $ n {A\ U A 2 ) is Poisson distributed with mean 

A(Ai U A 2 ) = A(^i) + A(4 2 ) - A(Ai n A 2 ). 

Because e~ A( - A ^ equals the retention probability h(ri,Wi) of grain (xi, ri,Wi) (see 
Proposition 3.1), we see that 

h 2 =h(r 1 ,w 1 )h(r 2 ,w 2 )e A ^ nA ^. 

The claim now follows after noting that 

A(AiHA 2 ) = A / \B ri+r (xi) ^ B r2+r (x 2 )\G r [w 1 V w 2 ,oo) F{dr). 

□ 

A key quantity for analyzing the covariance function of the thinned grain cover 
in Section 4 is the following function, which we shall call the retention covariance 
function. It is defined by 

q{u,r 1 ,r 2 ) = h 2 (u,n,r 2 ) - h(ri)h(r 2 ), (3.11) 

where h(r) denotes the weight-averaged retention probability defined in (3.3), and 
h 2 (u,ri,r 2 ) is the weight-averaged pair retention probability defined in (3.9). 

Lemma 3.3. The retention covariance function satisfies 

\q{u,r 1: r 2 )\ < h(r{) A h(r 2 ) 

for all u, ri,r 2 > 0. 

Proof. Fix a pair of reference grains (x\, r%, uui) and (x 2 , r 2 , w 2 ) having their centers 
at a distance u = \x\ — x 2 \ apart. Define a weight-dependent version of q by 

q(u,r 1 ,wi,r 2 ,w 2 ) = h 2 (u,rt, wt,r 2 ,w 2 ) - h(r 1 ,wi)h{r 2 ,w 2 ). 

We will first show that 

\q(u,r 1 ,w 1 ,r 2 , w 2 )\ < h{rx,wi), (3.12) 
by separately considering the following two cases: 

(i) If u < r\ + r 2 , then h 2 (u,ri,wi,r 2 ,w 2 ) is zero because the reference grains 
overlap, and (3.12) follows immediately. 
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(ii) If it > 7*1 + 7*2, then by borrowing the notation from the proof of Proposi- 
tion 3.2, we have by (3.10) that 

< h 2 (u,r 1 ,w 1 ,r 2 ,w 2 ) = P($(A 1 UA 2 ) = 0) < P(*(A t ) = 0) = h( ri , Wl ). 

As a consequence, 

-h(r 1 ,wi)h(r 2 ,w 2 ) < q(u,ri,wi,r 2 ,w 2 ) < h(ri, u>i)(l - h(r 2 ,w 2 )), 

from which (3.12) again follows. 

After integrating both sides of (3.12) over the weights, we see that \q(u, 7*i, r 2 )\ < 
h{r\). By symmetry, the same inequality holds with 7*1 replaced by 7*2, which proves 
the claim. □ 

3.5 Covariance function of the thinned grain cover 

Let us now consider the covariance function 



of the thinned grain cover X . 

Proposition 3.4. The covariance function of the thinned grain cover is given by 



where h is the weight-averaged retention probability defined by (3.3) and q is the 
retention covariance function defined by (3.11). 

Proof. Let us express the covariance function as 



where S±(z) is the probability that a single grain in $ simultaneously covers o and 
z, S 2 (z) is the probability that o and z are covered by distinct grains in < 3> th , and 
the volume fraction p t h can be viewed as the probability that an arbitrary reference 
point in IR rf is covered by some grain in <J? th . 

To write down an analytical expression for S\(z), recall first that by the hard- 
core property, the indicator function of X th can be written as 



k th (z) = P(o G X , z £ X th ) - P(o G X th ) P(z G X th ) 





B ri (o) n B r2 (x)\ q(\x - z\,r 1 ,r 2 )dxF(dr 1 )F(dr 2 ) 



(3.13) 



k th (z) = Si{z) + S 2 (z) - pth, 





(x,r,io)G<*> 



12 



where f y (Q; x,r,w) = l(y G B r (x))lj'^(x,r,w) is the indicator for the event that 
a grain (x,r,w) covers y and is contained in 3> th . Then 

Using Mecke's formula [SW08, Thm. 3.2.5] it's easy to see that 
S^z) = \ f \B r (p) n B r (z)\h(r) F(dr), 



where h(r) = J h(r,w) G r (dw). 

The probability that o and z are covered by distinct grains in can analo- 
gously be written as 

S 2 (z) = E ^ fo($;xi,r 1 ,wi)f z ($;x2,r2,w 2 ). 

(xi,ri,u>i)G<I> 
(x 2 ,r 2 ,ui2)e* 

Using the Slivnyak-Mecke formula [SW08, Cor. 3.2.3], it's not hard to verify that 
S 2 (z) = X 2 fff \B ri (o)nB r2 (x)\h2(\x-z\,r 1 ,r 2 )dxF(dr 1 )F(dr2), 



where h% is the pair retention probability defined by (3.8). The validity of the claim 
now follows after representing p t ^ using (3.7) and the identity |S ri (o)||B r2 (o)| = 
J \B ri {6) n B r2 (x) \ dx to note that 

Pth = A 2 Jff \B ri (o) n B r2 (x)\ h(r 1 )h(r 2 )dxF(dr 1 )F(dr 2 ). 

□ 



3.6 Two-point correlation function of thinned grain centers 

The two-point correlation function £th( z ) °f the random point configuration $* h = 
{x : (x, r, w) G $ th } of the thinned grain centers is defined as a function which 
satisfies 

Cov(**(A), $f(B)) = X 2 th J ^{x - y) dxdy 

for all disjoint and bounded measurable sets A, B C M d , assuming such function 
exists. This function, which in our case only depends on \z\, describes how much 
more (Cth( z ) > 0) or less (£th(z) < 0) likely it is to observe a point at a distance 
\z\ from a typical point, compared to observing a point in an arbitrary location. 
The two-point correlation function is related to the pair-correlation function g t ^ 
commonly used in statistics, via the formula t,th( z ) = 9th (z) ~ 1 ( e -g- [SKM95]). 
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Proposition 3.5. The two-point correlation function of the thinned grain centers 
is given by 

Cth(«) = 4-/ / q(\z\,r l7 r 2 )F(d ri )F(dr 2 ), (3.14) 

A th J®L+ 

where A t h is the thinned germ density defined in (3.5), and q is the retention co- 
variance function defined in (3.11). 

Proof. By using the Slivnyak-Mecke formula [SW08, Cor. 3.2.3] one can check that 
E<S>f(A)$f(B) = \ 2 I III h 2 (\x-y\,r 1 ,r 2 )F(dr 1 )F(dr 2 )dxdy. 

J A J B il+ JR+ 

for all bounded and disjoint A, B C M d , where h 2 is the weight-averaged pair 
retention probability defined in (3.9). On the other hand, Mecke's formula [SW08, 
Thm. 3.2.5] implies that 

E$f(A) = A / f h(r)F(dr)dx, 
J A JR+ 

where h(r) is the weight-averaged retention probability defined in (3.3). The claim 
follows by combining the above two formulas and recalling the definition of the 
retention covariance function (3.11). □ 

4 Long-range behavior of second-order statistics 

In this section we assume that the grain radius distribution F of the proposed 
Boolean model follows a power law with tail exponent a > d, by which we mean that 
the complementary cumulative distribution function F(r) = 1 — F(r) is regularly 
varying at infinity with exponent —a. In this case we can write 

F(r) = £(r)r' a , 

where the function t is slowly varying at infinity (see Appendix A for details). 
4.1 Asymptotic covariance 

The following result describes the covariance function of the thinned grain cover 
for thinnings where large grains have small retention probability. 

Proposition 4.1. Assume that the radius distribution F follows a power law with 
tail exponent a > d. Assume that the weight-averaged retention probability h{r) 
decays to zero as r — > oo, and that for any r\,r 2 > 0, the retention covariance 
function defined in (3.11) decays according to 

q(\z\,n,r 2 ) ~ q oc (r 1 ,r 2 )F(\z\)\z\ d as \z\ ->• oo. (4.1) 
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Then the covariance function of the thinned grain cover decays according to 

kth(z) ~ C-F(|z|)|,z| d as \z\ — > oo, 

where 

c = X 2 \B 1 \ 2 [ [ rf4q O0 (r 1 ,r 2 )F(dr 1 )F(dr 2 ). 

To prove Proposition 4.1 we need detailed results about the retention probabil- 
ities. The following lemma allows us to use dominated convergence on a part of the 
domain. 

Lemma 4.2. Assume that the radius distribution F follows a power law with tail 
exponent a > d. Then there exist constants c > and m > such that 

< q(\x - z\,n,r 2 ) < c\z\ d F{\z\) 

for all x, z £M. d and all ri, r 2 > such that \x\ < r\ + r 2 , \x — z\ > 1{r\ + r 2 ), and 
\z\ > m. 

Proof. Let c\ and u\ be the constants from Lemma B.2. Using the assumption 
that the function F follows a power law with tail exponent a, choose u 2 such that 
F(2/3r)/F(r) < 2(2/3)" a for all r > u x . Choose u 3 such that \ar d F(r) < 1 for 
all r > u 2 . Note that \x\ <r\+r 2 < ^\x — z\ implies \z\ < \x\ + \x — z\ < ||x — z\ 
and let m = max{ui, |u2, |^3}- Using Proposition 3.2 and the definition of m we 
have for all \z\ > m. 

q(\x - z\,ri,r 2 ) < exp J \B ri+r (o) n B T2+r (\x - z\)\F(dr)^j -1 

< exp ^Aci|x - z\ d F(\x - z\f 
<2Aci|x - z\ d F(\x - z\). 

Note that |x — z\ < \x\ + \z\ < ^\x — z\ + \z\ implies |x — z\ < 2\z\, and that F is a 
decreasing function. Now for c = 4(2/3) ~ a 2 d \ci and \z\ > m we have 

q{\x - z\,n,r 2 ) < 2Aci(2|z|) d F(||z|) < c\z\ d F{\z\). 



□ 



Lemma 4.3. Fix z £ M. d and define 

A(z) = {(x,rx,r 2 ) G R d x M + x R + : \x - z\ < 2(n +r 2 )}. 
Then the retention covariance function q satisfies 

\B ri (o) r\B r2 (x)\ \q(\x - z\,r 1 ,r 2 )\dxF{dr 1 )F(dr 2 ) 

A(z) 



<2\B 1 \ 2 ([ r d F(dr)] I [ r d F(dr) ] sup h(r), 

\Jr+ J \J\z\/6 / r>\z\/6 



\z\/6 ) r>\z\/6 

where h(r) is the weight-averaged retention probability defined by (3.3). 
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Proof. Define d[i as shorthand for dxF{dr\)F{dr 2 \ and denote the integrand by 
f z (x,ri,r 2 ). Observe that f z vanishes outside the set Aq = {(x,ri,r 2 ) : \x\ < 
i~l + r 2 }. Observe also that A(z) fl Aq C A\(z) U^C^), where Ai(z) = {(x,ri,r 2 ) : 
r,j > |z|/6}. As a consequence, 

fzdfi = / f z dfi < / f z dfi + / f z d/j, = 2 1 f z d^, 

A(z) JA(z)nA JAi(z) JA 2 (z) JA^z) 

where the last equality is due to the symmetry of f z with respect to its last two 
arguments. Recall that \q(\x — z\,ri,r 2 )\ < h(r±) by Lemma 3.3. Now 



A 1 (z) 



f z dfi< JJJ l^i/e^jCrOIB^CoJnBraCxJI/iCnJds^dnJFCdra) 
< J(z) sup h(r), 

r>\z\/6 



where 



J(z) 



l(|*|/6,oo)( r i) \B n (o) n B r2 (x)\dxF(dr 1 )F(dr 2 ) 



|Bi 



r d F(dr\ 



z|/6 



Proof of Proposition 4-1- By Proposition 3.4, we can write 

hh(z) 



where 



and where 



Ji(z) = [ \B r (o) n B r (z)\h{r)F{dr), 



H*) 
H*) 




fz(x,r 1 ,r 2 ) dxF(dri)F(dr 2 ), 




fz(x,r 1 ,r 2 ) dxF(dn)F(dr 2 ), 



denote the integrals of the function 



fz(x,ri,r 2 ) = \B ri (p) fl B r2 (x) 



q(\x - z\,ri,r 2 ) 



over the set 



F(\z\)\z\* 

A z = {(x, n,r 2 ) :\x-z\< 2(n + r 2 )} 



□ 
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and its complement, respectively. 

The integral I\{z) — > as \z\ — > oo by Lemma B.3, because h{r) — > as r — > oo 
by assumption. 

We will next show that I 2 (z) as well. We apply Lemma 4.3, to conclude 
that 

\h(z)\ <c 2 (\z\ d F(\z\))- 1 ( ^ r d F{dr) j sup h(r), 

\J\z\/6 J r>\z\/6 

where c 2 = 2\B±\ 2 f r d F(dr). The right side above tends to zero as \z\ — > oo, 
because h(r) — > as r — > oo, and because the integral on the right side above is 
asymptotically equivalent to constant multiple of \z\ F(\z\) by Lemma A.l. 

To analyze the limiting behavior of I%(z) as \z\ — > oo, note that assumption (4.1) 
and Lemma A. 2 imply that for any x,ri,r2, 

q(\x - z\,r 1 ,r 2 ) ~ q 00 {r 1 ,r 2 )\x - z\ d F(\x - z\) ~ qoo{ri,r 2 )\z\ d F(\z\). 

By the definition of A z , it thus follows that 

fz(x,r 1 ,r 2 )l A c(x,r 1 ,r 2 ) -> q oa (r 1 ,r2)\B ri (o) Pi Br 2 (a;)| 

as |z| — > oo. Moreover, by Lemma 4.2 there exists a constant C3 such that 

\fz(x,r 1 ,r 2 )lA%(x,r 1 ,r 2 )\ < c 3 \B ri (o) n -B r2 (x)| 

for all ar, r%,r 2 and all large enough z. Because the right side above is integrable 
with respect to dxF{dr\)F{dr 2 ), Lebesgue's dominated convergence theorem shows 
that 



lim h(z)= q 00 {r l ,r 2 )\B ri {o)^B r2 {x)\dxF{dr 1 )F{dr 2 ) 

|«|-K» JJJ 

= \Bi\ 2 [j qoo{r u r 2 )rir d F{dn)F{dr 2 ), 
which completes the proof of Proposition 4.1. □ 
4.2 Asymptotic two-point correlation 

Proposition 4.4. Assume that the radius distribution F follows a power law with 
tail exponent a > d and 

q(\z\,n,r 2 ) ~ g 00 (ri,r 2 )|z| d F(|z|) as \z\ -» 00. (4.2) 

Then 

£th(z) ~ c \ z \ F(\ z \) as \ z \ ~ * °°i 

where 



c=j 2 ~l j Qooin, r 2 )F(dr 1 )F(dr 2 ). 



th 
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Proof. Using (3.14) we can write 

fa(z) A 2 



where 



h(z)= // M^^Fidr^Fid^), 
J Ja z 

h(z)= ff fzi^r^Fidr^Fidrz), 



denote the integrals of 

, / \ q{\z\,n,r 2 ) 

over the set 

A z = {(n,r 2 ) E M+ x E+ : n + r 2 > \z\/2] 

and its complement, respectively. 

Observe that A z C A%(z) U A 2 (;z), where Ai{z) = {(ri 5 r 2 ) ■ > |^|/4}, and 
that |<jr| < 1 by Lemma 3.3. As a consequence, 

| g (|2|,ri,r2)|F(dri)F(dr 2 ) < (FxF)(A 1 (z)) + (FxF)(A 2 (z)) = 2F(|z|/4), 

which implies ii(#) — > as \z\ — > oo. 

Note that r 2 )lyic (n, r 2 ) — )■ (^(n,? - ^) by assumption (4.2) and the defi- 

nition of A z . By Lemma 4.2, f z (?i, r 2 )l J 4 c (?"i> ^2) is bounded for large z uniformly 
on n and r 2 . Lebesgue's dominated convergence theorem then shows that 



lim I 2 {z)= If q 00 {r l ,r 2 )F(dr)F(ds). 

z|-kx> 7 7 



□ 



5 Isolated grains retained 

In this section we study the thinning where only isolated grains are retained. In 
the general framework of Section 3, this is achieved by assigning unit weight to 
every grain, so that G r (dw) = 5\{dw). For nonrandom equally sized grains this 
corresponds to the classical Matern type I thinning. 

Theorem 5.1. Assume that the radius distribution F follows a power law with tail 
exponent a > d, so that 1 — F(r) = £(r)r~ a for some slowly varying function I. 
Then the thinned radius distribution is bounded by 

Fth(r) < Ae-Alftl^ 
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the covariance function of the thinned grain cover decays according to 

fcthO) ~ \c a , d pt h l(\z\)\z\-( a ~ d:} as \z\ -> oo, 
and the two-point correlation function of the thinned grain centers according to 

£ th O) ~ Xc a4 £(\z\)\z\- {a - d) as \z\ -> oo, 
where the constant c a ^ is given by (B.2). 

Proof. Because the weights are deterministic the retention probabilities have simple 
formulas 

fc(r) = /»(r,l)=exp^-A|.Bi| j (r + s) d F(ds)^j 



and 

h-2{\z\,ri,r 2 ) = h(ri)h(r 2 ) exp I A 



f \B ri+s (o)nB r2+s (z)\F(ds)) . 



The tail of the thinned radius distribution (3.6) is 

\ poo \ \ 

^th(r) = t — / h(s)F(ds) < —h(r) < — . 

\h Jr A th ^th 

To show the claim for the covariance and two-point correlation functions, we 
will use Proposition 4.1 and Proposition 4.4 respectively. For that we need to show 
that (4.1) holds. By Lemma B.l we have for the average intersection volume in h 2 
above 

\B n+s {o) n B r2+S (z)\ F{ds) ~ c\z\ d F{\z\ 



J 

Jr. 



Because the right hand side goes to zero as \z\ — > oo we can use the fact that 
lim^oo (e* — l)/t = 1 to obtain (4.1) with 

qoo(r 1 ,r 2 ) = h^ijh^Xc^d. 

Using Proposition 4.1 we find that kt\i{z) ~ c±\z\ ^(|^|). Using the formula 
for volume fraction (3.7) we also find the constant c\ = Ap^ h c a ^. Similarly by 
Proposition 4.4 we find that £th(-z) ~ C2|-z| d -F(|z|). With the help of germ density 
(3.5) we have c 2 = Xc a ^- □ 



6 Random grains retained 

Here we assume that each grain in the proposed Boolean model is assigned a random 
weight independently of the other grains, according to some continuous distribution 
function. The continuity ensures that there will be no tie breaks. Because the shape 
of the weight distribution does not affect the retention probabilities considered 
here, as long as it is continuous, we may without loss of generality assume that 
G r (dw) = l(o,i) (w)dw, the uniform distribution on (0, 1). Note that for nonrandom 
equally sized grains, this corresponds to the classical Materntype II thinning. 
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Theorem 6.1. Assume that the radius distribution F follows a power law with tail 
exponent a > d, so that 1 — F(r) = £(r)r~ a for some slowly varying function I. 
Then the thinned radius distribution decays according to 

^th(r) ~ (Athl^l)- 1 ^— £(r)r-(" +d ) as r -»■ oo, 
a + a 

the covariance function of the thinned grain cover according to 

k th (z) ~ didzDlzl-^-^ as \z\ ->• oo, 
and the two-point correlation function of the thinned grain centers according to 

£th(z) ~ c 2 e(\z\)\z\-( a - d ^ as |z| ->oo, 
for some c\,c 2 € (0, oo). 

Proof. The retention probability of a grain with radius r and weight w £ (0, 1) is 



h(r,w) = exp ^-A ^ oo)|# r+s (o)|F(<is)^ = exp(-A(l 



w)6(r)), 



where 



b(r)= [ \B r+s (o)\F(ds). 



The weight-averaged retention probability thus equals 
fc(r) = / e- x ^ w ^dw = 1 



A6(r) 

Because 6(r) — > oo as r — > oo, it follows that lim^oo h(r) = 0. The first condition 
of Proposition 4.1 is thus satisfied. 

Note that b(r) ~ |i? r | as r — > oo, which implies that 

h{r) ~ (Afo(r))- 1 ~ (AjSijr^)- 1 . 

By Lemma A. 4, the thinned radius distribution (3.6) is 

\ /*00 /'OO 

*th(r) = — / h(s)F(ds) ~ (AthlBil)- 1 / S ~ d F(cte). 

^th Jr Jr 

Furthermore by Lemma A.l 

^(^-(Athl^l)- 1 ^-^^). 

a + d 

The pair retention probability equals 

h2(.\z\,n,r2,wi,w 2 ) 
=h(ri,wi)h(r 2 ,w 2 )exp ( A / V tu 2 , oo)|B ri+s (o) n B r2+s (z)\F(ds) 

= exp ( - A(l - wi)6(n) - A(l - w 2 )b(r 2 ) + A(l - w% V iu 2 )a 2 (ri, r 2 )J , 
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where 



a z (ri,r 2 )= / \B ri+s (o) D B r2+S (z)\ F(ds). 



From this expression we see that the retention covariance function defined in (3.11) 
equals 



g(|4 ri,r 2 ) = f j e - A6 (n)(l-^i) e -A6(r 2 )(l-™ 2 ) ^(l-toiVteaMn.tT,) _ ^ 



7 1 

J 

As 1 2 1 oo, Lemma B.l shows that the term in parentheses above is asymptotically 
equivalent to 

e A(l-«iV« 3 )a.(n,r a ) _ j „ A(l - Wi V W 2 )c a>d F(\z\) \z\ d . 

With the help of the bound |e* — 1] < (e — l)t for i € [0, 1], we may use dominated 
convergence to conclude that 

q(\z\,ri,r 2 ) ~ q oc (r 1 ,r 2 )\z\ d F(\z\) as \z\ -> oo, 

where 



jo Jo 

Now by Proposition 4.1 it follows that 

^th(^) ~ ci-P(|^|) as \z\ — > oo, 

where 

ci=\ 2 \B 1 \ 2 [ [ rfrlq 00 (r 1 ,r 2 )F(dr 1 )F{dr 2 ). 

The constant c\ is finite because <Zoo(rL, ^2) < Ac QjC ; for all ri,r 2 . The fact that ci 
is strictly positive is easily seen by inspecting the expression of qoo{ri-,r 2 ). 
Similarly, Proposition 4.4 shows that 

&h(*) ~ c 2 \z\ d F(\z\), 

where 

C2 = TT / Qoo{ri,r 2 )F(dri)F(dr 2 ). 
\h JR+ J1+ 

The finiteness and strict positivity of c 2 follow by similar reasoning as for c\. □ 
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7 Large grains retained 

A thinning which favors large grains is obtained by letting the weight of each grain 
be equal to its radius, so that G r (dw) = S r (dw). 

Theorem 7.1. Assume that the radius distribution F follows a power law with tail 
exponent a > d, so that 1 — F(r) = £(r)r~ a for some slowly varying function I. 
Then the thinned radius distribution decays according to 

Fth(f) ~ t — £(r)r~ a as r — > oo, 
Ath 

the covariance function of the thinned grain cover according to 

k th (z) ~ \c a4 (l - p th ) 2 i{\z\)\z\- {a - d) as \z\ oo, 
and the two-point correlation function of the thinned grain centers according to 

6h(s) ~ Xc a , d £(\z\)\z\'^ as \z\ oo, 
where the constant c ai d is given by (B.2). 

Proof. Because the weight of each grain is equal to its radius, the weight- averaged 
retention probability h(r) is equal to h(r, w) with w taking on the value r. By 
Proposition 3.1, the retention probability is given by 

h(r) = exp^-A^ \B r+s \ l [f . )00 )(s) F(ds)J . 

Because the integrand above tends to zero as r — > oo, and the integrand is bounded 

by the i ? ((is)-integrable function \B2 S \, dominated convergence implies that lim,.-^ h(r) = 

1. By Lemma A. 4, the tail of the thinned radius distribution (3.6) satisfies 

Ah(r) = f - / h(s)F(ds) ~ / F(ds) = ±F(r). 

A th J r A th J r A th 

To analyze the long-range behavior of k t h(z) and £th( z )i ^ us ^ IS ^ investigate 
the long-range behavior of the retention covariance function q(\z\, r\, T2) defined 
by (3.11). Using Proposition 3.2, we find that 

q(\z\, ri ,r 2 ) = h(n)h(r 2 ) (l(n + r 2 < \z\)e T ^ r ^ - l) , (7.1) 

where 

r(\z\,r 1 ,r 2 ) = X / l[ riVr2)00 )(s)|S ri+a (o) fl B r2+S (z)\ F(ds). (7.2) 

When \z\ > 3(ri + r 2 ), we may replace the region of integration above with the full 
positive real line, so that with the help of Lemma B.l we find that 

T(\z\, n ,r 2 ) = A / \B ri+s (o)DB r2+s (z)\F(ds) ~ \c a4 \z\ d F(\z\), 



22 



as \z\ — > oo. Because & — 1 ~ t for small t, we conclude using (7.1) that 

q(\z\, ri ,r 2 ) ~ (?00 (r 1 ,r 2 )|z[ d F(|z|), (7.3) 

where 

9oo(ri,r 2 ) = \c ayd h(ri)h(r 2 ). (7.4) 

The claim for the two-point correlation function £th(-z) now follows by using Propo- 
sition 4.4, after noting that the constant in Proposition 4.4 is 



qoo(ri,r 2 )F(dn)F(dr 2 ) = \\7*c a4 ( A / h(r)F(dr) ) = \c a4 . 



We will now move on to the part concerning the covariance function fcth(-z) °f 
the thinned grain cover. Note that because h{r) does not vanish as r — > oo, we 
cannot use Proposition 4.1 to deduce the long-range behavior of kth(z)- Instead, 
we will proceed by directly analyzing the integral building blocks of feth(^) i n high 
precision. Let us start by rewriting (3.13) as 



kth(z) = AIo(^) + \ 2 {h{z) + h{z) + h(z)), 



where 



h(z)= / \B r (o)nB r (z)\h(r)F(dr), 
Jr + 

I j (z)= III \B ri (o)nB r2 (x)\q(\x-z\,r u r 2 )dxF(drx)F(dr 2 ), j = 1,2,3, 




A: 



and 



A\ = {(x, ri, r 2 ) : \x\ < n + r 2 , r% + r 2 < \x - z\/2}, 

A 2 = {(x, ri,r 2 ) : \x\ < r\ + r 2 , \x - z\/2 < r\ + r 2 < \x - z\}, 

A\ = {(x, ri,r 2 ) : \x\ < r\ + r 2 , \x — z\ < r\ + r 2 }. 

The first term Iq{z) < c a ^\z\ d F(\z\) by Lemma B.l. Note that the integrand 
in Iq(z) vanishes for r < \z\/2 so that 

I (z) > inf h(r) f \B r (o)nB r (z)\F(dr). 

r>\z\/2 J R+ 

Using Lemma B.l and that h{r) — > 1 as r — > oo we conclude that 

h{z) ~ c a4 \z\ d F(\z\). (7.5) 

Next, we will prove that 

h{z)~\- l c a4 pl h \z\ d F(\z\). (7.6) 
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By Lemma 4.2, the function (x,ri,r 2 ) h-» |fp§fefr^ 0^ r i ? r 2) is positive and 
bounded by a constant which does not depend on z. Because \B ri (o) n 5 r2 (x)| 
is integrable with respect to dxF(dri)F(dr 2 ), Lebesgue's dominated convergence 
theorem shows that 



l ™jM^= III l^»nB r2 (x)| ( lim ^-J^ Uf ) dxF^F^). 



Using (7.3) and the definition of Af, the limit on the right equals g o('*i) ^2)- Plug- 
ging in the expression (7.4) for qoo( r i, r 2) an d recalling the formula (3.7) for the 
volume fraction of the thinned grain cover p t ^, we find that 

lim i F\ = Xc ^ d III \B r MnB r2 (x)\h{rx)h{r 2 )dxF(drx)F(dr 2 ) 

|z|^oo \z\ d F(\z\) J J J 

= Ac, 



, d J J iBrfrfrthinM^Fidr^Fidrz) 



= A 1 c a4 p 2 th , 

which proves the validity of (7.6). 
Now we will prove that 



- as |*| -> 00. (7.7) 



,\*F(\z\) 



First, using the bound \B ri+s (o) n \B T2+S {x — z)\ < \B ri+s \ < \B2 S \ for s > r\ V r 2 , 
we find that the function r defined in (7.2) is bounded by 



t{\x - z\,n 

Observe next that 



,r 2 )<X\B 1 \2 d J l [riVr2t0o) (s)s d F(ds). 



\z\ < 3(n + r 2 ) < 6(ri V r 2 ) 

for all (x,ri,r2) G vlf, so that n V r 2 is large when |z| is large. As a consequence, 
we see by Lemma A.l and Lemma A. 3 that for all (x,ri,r 2 ) G A\ and all large 
enough z, 

t(\x - z\,r u r 2 ) < 2A| J Bi|2 d (r 1 \J r 2 ) d F{n Vr 2 ) 
<4A|Bi|2 d (|z|/6) d #(|*|/6). 

Because e* — 1 < (e — l)i for i G [0,1], formula (7.1) combined with the above 
inequality shows that for all (x,ri,r 2 ) G A 2 and all large enough z, 

< q(\x - z\, ri ,r 2 ) < Cl \z\ d F(\z\/6), 

where c\ = 4(e — l)X\Bi\3~ d . Therefore, 

< h{z) < d\z\ d F(\z\/6) /// \B ri (o) r\B T2 {x)\dxF{dr 1 )F{dr 2 ). 
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Note that A% C A\\ U A z 22 where A z 2i = {(x,n,r 2 ) : r» > |z|/6}, % = 1,2. By 
symmetry of the integrand with respect to n and r2, 



< I 2 {z) < 2 Cl \z\ d F(\z\/6) I I I \B ri (p) n B r2 (x)\ dxF{dr l )F{dr 2 ) 

1*1/6 ■ 



2ciH d J F(|z|/6)|5 1 | 2 r d F{dr)\ / r d F(d 



OC 



1*1/6 



which shows the validity of (7.7). 
It remains to be shown that 

I 3 (z)~- 2 ^c atd \z\ d F(\z\). (7.8) 

To do that, we first fix a small e 6 (0, 1/4). Note that by formula (7.1), the retention 
covariance function equals q(\x — z\, ri, r 2 ) = —h(ri)h(r 2 ) for (z, 7i,r 2 ) G Note 
also that, for fixed ri and r 2 the x-slice of A| is 

{x : (x,ri,r 2 ) G A 3 } = £ ri+r2 (2:) nB n+f2 (o). 

Because |jB n (o) nB r2 (a;)| vanishes for x outside B ri+T2 (o), we may represent Is(z) 
according to 




h{z) = ~ /// \B ri (o)nB r2 (x)\h(r 1 )h(r 2 )dxF(dr 1 )F(dr 2 ) 

\B n (o) n B r2 (x) | dx h(r x )h{r 2 ) F(dr 1 )F(dr 2 ) 




C| J B ri + r2 (z) 



where 

Cl = {(n,r 2 ):\z\<2( ri + r 2 )}. 
Next we split ^(z) into three parts 

h{z) = -(hi(z) + I 32 (z) + / 33 (z)), 

where 



7 3iW= // / |S ri (o)n J B r2 (x)| ( ix/ i (r 1 )/ i (r 2 )F(dr 1 )F(dr 2 ) 

JB ri+r2 (z) 



for j = 1, 2, 3 and 



Cfi — c 311 U C 312 

Cfn = {(n,r 2 ) : < ri < e\z\, \z\/2 < r 2 } 

C312 = {(n,r 2 ) : < r 2 < e|z|, |z|/2 < n} 
C32 = C331 u C332 

C 3 Z 2l = {(ri,r 2 ) : < n < e|^|, M/2 - n < r 2 < \z\/2} 

Cf 22 = {(n,r 2 ) : < r 2 < e\z\, \z\/2 - r 2 < n < \z\/2} 
CI 3 = [e\z\,oo) 2 nC$. 
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A change of variables shows that 

/ \B Tl (o) D B r2 (x)\ dx = \B r , 2 (z) n B ri J rr2 (x) \ dx, 

J B ri+ r 2 (z) J B ri (o) 

so that we can express the integral -^(z) more conveniently as 



hi{z)= // / \B r2 {z)^B ri ^ 2 {x)\dxh{r^h{r 2 )F{dr 1 )F{dr 2 ). 

J JC^ JB ri (o) 

By symmetry we can write -^31(2) = 2/311 (z), where 1311(2;) is a modification of 
131(2) with the region of integration CJ^ replaced by C| n . To analyze the long- 
range behavior of 1311(2), let us split it according to 1311(2) = J\{z) + J 2 (z), where 

J 1 (z) = / / / \B r2 {z)nB r2 {o)\dxh{r2)F{dr 2 )h{r 1 )F{dr l ), 

JO J\z\/2 J B ri (o) 

and where J 2 (z) = ^311 (z) — J\{z). Because the integrand of Ji(2) does not depend 
on x, we can rewrite the integral as 

J x [ z ) = ^jf ' ' \B r \ h(r)F(dr)j ^ J ^ \B r (z) n B r (o)\h(r)F(dr)j . 

The first integral on the right satisfies 

/' \B x \r d h(r)F{dr) ~ f |S 1 |r d /i(r)F((ir) = A _1 pth, 
Jo Jr+ 

where p t h is the volume fraction of the thinned grain cover given by (3.7). Notice 
that, because the intersection in the second integral vanishes for r < \z\/2, we can 
apply (7.5) to conclude that 

Ji(z) ~ \- l p^c a4 \z\ d F{\z\). 

The rest of the proof constitutes of showing that the remaining three parts of 
13(2) are negligible. We start by showing that J 2 > and 

- |Bl ' 2 ( (1 +2t)d ~ ') (/ rdF(dr) ) 

First we need a bound for the difference of the intersections in J 2 (2). Fix x £ B ri (o) 
and (ri,r 2 ) € C| 11 . Because \x\ < r\, we have B r2 {o) C B ri+r2 (x), which implies 
that the integrand in J 2 (z) is bounded by 



< 



< 



< 



B r2 (z) n B ri+r2 (x) I - \B r2 (2) n B r2 (o) 
B r2 (z)n(B ri+r2 (x)\B r2 (o))\ 
B, n+r2 (x) \ - \B r2 (o)\ 
B l \{{l+r 1 /r 2 ) d -l)r d 2 
B l \((l + 2e) d -l)r d , 
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where the last inequality is due to n < e\z\ and \z\/2 < r 2 . This bound and 
h(r) < 1 now imply that 



0<J 2 {z)< [ [ [ [S 1 [((l + 2e) d -l)r^x/i(r 2 )F(dr 2 )/ i (r 1 )F( ( ir 1 ) 

JO J\z\/2 J B ri (o) 

< |5i| 2 ((l + 2e) d - 1) r d F(dr)^ r d F(dr)^j . 



Now using Lemma A.l proves the claim. 
We will now show that ^32 > and 



^W$k s 2 (/ /F(dr) ) |Bl ' 2 [ (1 - 2£) "° - 1 



" -2"-". 



a — d 



By symmetry 132(2) = 2/321(2), where 

1321(2)= // / iB^WnB^t^ldx^^^JF^OF^), 

^321 ^v-^ ~| _ ^'2 (^0 

Note that C| 21 C [0, e[z|] x [|z|/2(l — 2e), \z\/2\. Also approximating i? ri+r2 (z) by 
M. d and recalling that /i(r) < 1 we have 



f \z\e r \z\/2 

Imi{z) < / / \B 1 \ 2 rir d F(dr 2 )F(dr 1 ) 

JO J(\z\/2)(l-2e) 



•\z\e r\z\/2 

'(|*|/2)(1-2,1 

nOO 



r d F(dr) - / r d F(dr) 

(|*|/2)(l-2 e ) J\z\/2 



Now using Lemma A.l implies the claim. 

For the last part 133(2) we have first a simple bound 



733(2) < / / / |B n (o)n5 ra (a;)|dar/i(ri)/»(r2)F(dri)F(dr 2 ) 

\z\ J B ri+r2 (z) 



£\Z\ J € 
OO /"OO 



< 



/*oo /*oo 

/ / l^il 2 r d r d F(dn)F(dr 2 ) 

J e\z\ J e\z\ 

/ 00 \ 2 

y 6 [*i / 



Using Lemma A.l this bound implies that 

^33(2) 

, 1 . — j r — > as 2 — > 00. 
\z\ d F(z) 1 1 

Adding together all the parts of ^3(2) we have 

,. h{z) Pth 
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and 

|z|^oo |z| a F(|z|) A 

where 



r d F{dr) ) -^— [ 2a ~ d { { (1 + 2 e) d - l) + ( C 1 - 2e) d ~ n 



,5(e) = 2|i? 1 | 2 (7<, 

Letting e — > shows the validity of (7.8) and concludes the proof. □ 



8 Small grains retained 

In this section we study a thinning which favors small grains. This thinning is 
obtained by setting the weight of each grain to the inverse of its radius, so that 
G r {dw) = 5\/ r {dw). The following theorem shows that the thinned radius distri- 
bution and the key second-order characteristics decay rapidly to zero, regardless of 
the tail behavior of the original radius distribution F. Note that here, unlike in 
Theorems 5.1-7.1, there is no need to assume anything on the shape of the radius 
distribution F. 

Theorem 8.1. Assume that the radius distribution F satisfies f r d F(dr) < oo. 
Then the thinned radius distribution is bounded by 

F th (r) < ±e~WK 
the covariance function of the thinned grain cover by 

and the two-point correlation function of the thinned grain centers by 

|&h(*)|<e-WI*l d 
for all large values of r and z, where c G (0, oo). 

Proof. Now the weight-averaged retention probability h(r) is equal to the retention 
probability of a reference grain with radius r and weight 1/r. Using Proposition 3.1 
we find that 



h(r) = exp f-Ajf \B r+s (o)\F(ds)^j 



From this expression we see that h decreases monotonically to zero as r grows, and 
that h(r) < exp(— ^X\Bi\r d ) for all large enough r so that F(r) > 1/2. Proposi- 
tion 3.2 further shows that the weight-averaged pair retention probability equals 



h2(\z\,r 1 ,r 2 ) = h(n)h(r2) exp I A / l [0jriAr2 ](r)|.B n+r (o) n B r2+r (z)\F(dr) 
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for \z\ > T\ + T2- 

To analyze the covariance function of the thinned grain cover, recall that 



k th (z) = X I \B r (p) n B r {z)\h(r)F(dr) 

+ A 2 Iff \B ri (o)nB r2 (x)\q{\x-z\,r 1 ,r 2 )dxF{dr l )F(dr2), (8.2) 



where q(u,n,r2) = h2(u,rx,r2}—h(ri)h(r2)- Because \B r (o)nB r (z)\ < \Bi\r d l(r > 
|z|/2), the first term on the right side of (8.2) is bounded from above by 



A|5i| ( f r d F(dr)\ sup h(r) 

\ J\z\/2 I r>\z\/2 



z\/2 J r>\z\/2 

Note that q(\x — z\,r\,T2) vanishes for \x — z\ > 2{r\ + r 2 ), because the integral 
in (8.1) vanishes for \z\ > 2(r% + r2). This is why the integration in second term 
in (8.2) can be restricted to the set A(z) = {(x, r±,r2) '■ \x — z\ < 2(ri + ^2)}. Now 
using Lemma 4.3, the absolute value of second term in (8.2) is bounded from above 
by 



) J [ / r d F(dr) ) sup h{r) 
J \J\z\/6 J r>\z\/6 



2X 2 \B 1 \ 2 / r d F{dr) 

As consequence, we find that 

|feth(- 2 )| < (Ami + 2A 2 mf) sup h(r), 

r>\z\/6 

where m\ = \Bx\j r d F(dr) is the mean volume of a grain. Therefore, 

for all large enough z such that i^(|^|) > 1/2. A similar analysis can be carried out 
for the two-point correlation function. 

Typical radius has tail probabilities (3.6) 

\ poo \ \ 

F th (r) = -f- / h(s)F(ds) < -— h{r)F(r) < -^- e -^ Bl \ rd 

Ath Jr A t h A t h 

for all large enough r so that F(r) > 1/2. □ 



9 Conclusions and future work 

Boolean models consisting of randomly sized spheres in M. d are long-range dependent 
if the sphere radii follow a power-law distribution with tail exponent a € (d, 2d). 
We studied second-order statistical properties of four hard-core germ-grain models 
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which are obtained from such Boolean models using a natural weight-based thin- 
ning mechanism. We found that a thinning which favors large grains preserves the 
power-law covariance decay of the proposed Boolean model, whereas a thinning 
which favors small grains does not. The power-law covariance decay is also pre- 
served under the thinning where only isolated grains are retained (Matern type I), 
and the thinning where retention is determined by independent weights (Matern 
type II). The germ-grain model obtained by the Matern type I thinning is an in- 
teresting example of a homogeneous hard-sphere model where typical spheres have 
exponentially small sizes but the covariance function decays slowly according to a 
power law. 

To keep the notation simple and the paper easy to read, the analysis carried 
out in this article was restricted to spherical grains. However, we believe that 
this assumption can be easily relaxed to some extent following standard techniques 
of stochastic geometry. Another interesting open problem is to investigate how 
thinnings affect covariance decay properties in the light-tailed setting where the 
grain size distribution in the proposed Boolean model is assumed to decay rapidly. 



A Regular variation 



A measurable function / : R + — > R is called regularly varying with exponent 7 € R 
if it is positive for all large enough input values and for all a > 0, 

/(<) 

as t — > 00. A regularly varying function with exponent zero is called slowly varying. 
For a good overview on the theory of regular variation, see for example [BGT87]. 
In this section we will summarize some key properties of regularly varying functions 
which are needed in the analysis. The first one is a consequence of the Karamata's 
theorem [BGT87]. 

Lemma A.l. Let F be distribution function on R + with a regularly varying tail of 
exponent a > p. Then for any constant a > 0, 



r p F(dr) 



a 



a — p 



-(«-p) _p| 



XjX 1 ' 



as x — >• 00. 



Lemma A. 2. Assume that i is slowly varying. Then for any zq £ 

£(\z-z \) ~£{\z\) as\z\ 



00. 



Proof. Fix zq £ R d , and write \z — zq\ = a z \z\, where a z = \z — zq\/\z\. Because 
a z — > 1 as \z\ — > 00, we can fix m such that a z £ [1/2,3/2] for \z\ > m. Now, for 
any z £ R d such that | jar | > m, it follows that 



\Z ~ Zq\ 



£(a z \z\) 



< sup 

06(1/2,3/2] 



i(a\z\) 
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The right side above tends to zero as \z\ — > oo because €(o|^|)/-£(|^|) — > 1 locally 
uniformly with respect to a [BGT87, Thm. 1.2.1]. □ 

Lemma A. 3. Assume that f is regularly varying with exponent —7 < 0. Then 
there exists a constant u > such that f(y) < 2/(x) for all y > x > u. 

Proof. By the uniform convergence of regularly varying functions [BGT87, Thm. 1.5.2], 
f{\x)/f{x) — > A~ 7 uniformly for A > 1 as x — > 00. This implies that we can find 
u > such that /(Ax) < 2/(x) for all x > u and all A > 1. Now because y > x, we 
have 



f(y) < sup /(Ax) <2/(x). 
A>1 



□ 



Lemma A. 4. Let F a probability measure on K + , and let f and g be bounded 
positive functions on M + such that f(r) ~ g(r) as r — > 00. Then, 



f(s)F(ds) 



g(s)F(ds) as r 



00. 



Proof. 



rmF(ds 



J r °° g(s)F(ds) 



r (1 - SO s(s) F (ds) 



L°° g(s)F(ds) 



< sup 

r<s 



□ 



B Intersections of distant balls 

Lemma B.l. Let F be a probability distribution on M+ which follows a power law 
with tail exponent a > d. Then for any r\,ri > 0, 

f \B ri+r (o)DB r2+r (z)\F(dr) ~ c a4 F{\z\)\z\ d as \z\ -> 00, (B.l) 

where 

c a , d = / \B r (o)nB r (e 1 )\ar- a - 1 dr, (B.2) 
Jo 

and e\ is the first unit vector in the standard basis ofM. d . 

Proof. Because the Lebesgue measure is rotation-invariant, we may assume without 
loss of generality that z = ue% for u > 0. Fix ri,T2 > 0, and denote left side of 
(B.l) by I(u). We will prove the claim by first showing that 

Ji(«) = / \B r (o) n B r {ue 1 )\F{dr) ~ cu d F(u), (B.3) 
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and then showing that the remainder / 2 (w) = I(u) — I\(u) tends to zero faster than 
u d F{u) as u — > oo. 

To prove (B.3), let F u be the distribution of a random variable obtained by 
dividing a -F-distributed random variable by u, so that F u [r) = F(ur). Then a 
change of variables shows that 



h{u) = u d \ \B r/u {o) n B r/u ( ei )\F(dr) = u a / <f>(r) F u (dr), 

where (j>(r) = [-B r (0) Pi B r (e\)\. Because <fi is continuous, r~ d cj)(r) < \B\\ for all 
r > 0, and (j){r) = for r < 1/2, we may apply [KLNS07, Lemma 2] (with p = d, 
7 = a, q = a + 1) to obtain 

/ <j>(r) F u (dr) ~ F(u) [ </>(r) ar'^dr = cF(u), 
Jr + Jo 

which implies the validity of (B.3). 

To show that l2(u) tends to zero faster than u d F(u), note first that for all 
u > 2(n + r 2 ), 

/•oo / \ roo 

h{u)= / r d ip u (r)F(dr) < sup ^ u (r) / r d F(dr), (B.4) 
Ju/4 Vr>w/4 I Ju/4 

where 

ip u (r) = \B 7 , l/r+1 (o) n B r2 / r+1 (e 1 u/r)\ - \B 1 (o) D Bi(eiu/r)\. 

The equality in (B.4) follows because ip u { r ) = when u > 2(ri + r 2 ) and r < u/4. 

Note that by Lemma A.l, the integral on the right side of (B.4) is asymptotically 
equivalent to 

' ^ r d F(dr) ~ f — ^— J 4 a ~ d F(u)n d . 



w/4 \a-d^ 
In light of (B.4), it hence suffices to show that 

sup V«( r ) (B.5) 

r>u/A 

as u — > oo. This will be done by inspecting the geometry of ip u . Because the 
intersection of the unit-balls above is a subset of the intersection of the larger balls, 
we can bound the nonnegative function ifj u using the annuli around the unit-balls, 
so that 

Mr) = \{B ri/r+1 (o) H B r2/r+1 (e lU /r)) \ {B x (o) n B x (e x u/r)) \ 



< 



\ B n/r+l(o) \ Bx(o)\ + \B r , 2/r+1 (eiu/r) \ B^au/r)] 
\B x \{{r x /r + l) d - 1) + [Ba|((r 2 /r + l) d - 1) 



<2|5 1 |((^ + 1) -1 



Because this bound is valid for all n, we conclude (B.5), and the proof is complete. 

□ 
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Lemma B.2. Let F be a probability distribution on M+ which follows a power law 
with tail exponent a > d. Then there exist constants u > and c > such that 

/ \B n+r (o) nB r2+r (z)\F(dr) < cF(\z\)\z\ d (B.6) 
Jr + 

whenever \z\ > u and r\ + r 2 < \z\/2. 
Proof. Observe first using Lemma A.l that 

r d F{dr) ~ Cl \z\ d F(\z\), 

z\/4 

where c\ = A a ~ d a/{a — d). Hence, we may fix a constant u > such that 

POO 

/ r d F(dr) < 2ci\z\ d F(\z\) (B.7) 
7 1*|/4 

whenever \z\ > u. 

Assume now that \z\ > u and r\ + r<i < \z\/2. In this case the intersection on 
the left side of (B.6) is nonempty only when r > |z|/4. For any such r > \z\/A, a 
crude estimate shows that 

|5 n+r (o)nB r2+r (^)| < \B ri+r (o)\ < \B 1 \(r 1 +r 2 + r) d , 

which together with the inequality r\ + r 2 < \z\/2 < 2r shows that 

\B n+r (o)nB r2+r (z)\ < 3 d \B 1 \r d . 

As a consequence, 

r roc 

/ \B ri+r (o) n B r2+r (z)\F(dr) < 3 d \B 1 \ / r d F(dr), 

Jr + ' 

so that by virtue of (B.7), the claim holds for c = 2ci3 d |-Bi|. □ 

Lemma B.3. Let F be a probability distribution on which follows a power law 
with tail exponent a > d, and let h be an arbitrary positive function. Then there 
exist constants u > and c > such that 

(IzfFdzl))- 1 [ \B r (o)DB r (z)\h(r)F(dr) <c sup h(r) 

JR + r>\z\/2 

for \z\ > u. 

Proof. By using Lemma B.2, fixing the constants u and c as in the lemma, and 
noting that the integrand vanishes for r < \z\/2, we see that 



\B r (o) n B r (z)\h(r)F(dr) < h*(z) / \B r {o) n B r (z)\F(dr) < ch*(z)F(\z\)\z\ d 

Jr+ 

for all \z\ > u, where h*{z) = sup J ,>| z |/ 2 h{r). □ 
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